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Abstract 



(N 

■^j- We study a general class of supersymmetric AdS4 x Y-j solutions of M-theory that 

have large N dual descriptions as J\f = 2 Chern-Simons- matter theories on S 3 . 
The Hamiltonian function Km for the M-theory circle, with respect to a certain 
contact structure on Yj, plays an important role in the duality. We show that an 
M2-brane wrapping the M-theory circle, giving a fundamental string in AdS,4, is 
supersymmetric precisely at the critical points of Km, and moreover the value 
of this function at the critical point determines the M2-brane action. Such a 
configuration determines the holographic dual of a BPS Wilson loop for a Hopf 

S^ circle in S s , and leads to an effective method for computing the Wilson loop on 

both sides of the correspondence in large classes of examples. We find agreement 
in all cases, including for several infinite families, and moreover we find that the 
image Um^Xi) determines the range of support of the eigenvalues in the dual 
large N matrix model, with the critical points of Km mapping to points where 
the derivative of the eigenvalue density is discontinuous. 
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1 Introduction and summary 

Over the last few years our understanding of the AdS4/CFT 3 correspondence, par- 
ticularly in M-theory, has improved considerably. Broadly speaking, this has involved 
developments on two fronts. Firstly, we now have large classes of very explicit examples 
of dual pairs; that is, gravity backgrounds for which we have some precise description 
of the dual superconformal field theories. Secondly, there are new quantitative tests of 
these conjectured dualities, based on supersymmetric localization in the field theories. 
The aim of this article is to extend this quantitative analysis further, by examining the 
computation of certain BPS Wilson loops on both sides of the correspondence. In the 
process we will also understand how other structures are related via the duality. 



Starting with the seminal work of [T] we now have large classes of supersymmetric 
A0IS4 x Y 7 gravity backgrounds of M-theory that are associated with particular (2 + 1)- 
dimensional supersymmetric gauge theories, typically Chern-Simons theories coupled to 
matter, that are believed to have a dual superconformal fixed point. The construction 
of the UV gauge theory usually relies on a dual description in terms of type IIA string 
theory, which in turn involves a choice of M-theory circle U(1)m acting on Y?; different 
choices can lead to different UV gauge theories that flow to the same IR superconformal 
fixed point. In [Tj the highly supersymmetric case where Y 7 = S 7 /Z&, equipped with 
its round Einstein metric and with N units of flux through this internal space, was 
related to a large N dual description as an M = 6 superconformal U(N) x U(N) 
Chern-Simons- matter theory (the ABJM theory), with fceZ being the Chern-Simons 
coupling. Here Z fc C U{1)m, with the M-theory circle action by U(1)m being the Hopf 
action on S 7 , so that S 7 /U(1)m — CP 3 . There are now many families of examples of 
a similar type [2] -[IE], generally with J\f > 2 supersymmetry, in which Y? is a Sasaki- 
Einstein seven-manifold and the dual description typically involves supersymmetric 
Chern-Simons-matter theories whose gauge groups are products of unitary groups, and 
with matter in various representations (bifundamental, fundamental, adjoint). There 
are also examples in which AdS4 x Y 7 is a warped product, with non-trivial four-form 
flux on non-Einstein Y 7 (obtained thus far either by marginal [T9] or relevant [2U] [25] 
deformations of Einstein examples). 

Quantitative tests of these conjectured dualities arise by putting the (Euclidean) 
field theories on a compact three-manifold. The simplest case, in which this three- 
manifold is taken to be S 3 equipped with its round metric, was studied in J2JH 1271 128] . 
This can be done for a completely general M = 2 supersymmetric gauge theory, in 
such a way to preserve supersymmetry. Moreover, using a standard argument [29] one 
can show that the path integral, with any BPS operator inserted, reduces exactly to 
a finite-dimensional matrix integral. This implies that the VEVs of BPS operators 
may be computed exactly using a matrix model description, with the large N limit of 
this then expected to reproduce certain supergravity results. In practice this has been 
used to compute the free energy F (minus the logarithm of the partition function) 
on both sides of the correspondence [30] [37], where on the supergravity side this is 
proportional to iV 3//2 with a coefficient depending only on the volume of Y 7 ]q 

It is natural to try to extend these results further, by inserting non-trivial BPS 
operators into the path integral, computing the corresponding large N behaviour in 
the matrix model, and comparing to an appropriate dual semi-classical supergravity 
computation. In the original papers on the ABJM theory [261 EDI ESI ESI EH SI] 



-"^For a general AdS4 x Yj solution this is the contact volume of I7, rather than the Riemannian 
volume, as we shall review in section [3] 



the supersymmetric Wilson loop for the gauge field around a Hopf circle S 1 C S 3 was 
studied. This is 1/2 BPS, and is readily computed in the large N matrix model [261 130]. 
Generally speaking, one expects such a Wilson loop to be dual to a fundamental string 
when viewed from a type IIA perspective [42], with the Euclidean string worldsheet 
having boundary on the Hopf S 1 at conformal infinity. Semi-classically, more precisely 
this will be a supersymmetric minimal surface £2 m Euclidean AdS4, with the VEV 
calculated via the regularized area of the string worldsheet. Such a string must then 
be pointlike in the internal space, and for the ABJM theory this is CP 3 = S 7 /U(1)m- 
Equivalently, this IIA string lifts to an M2-brane wrapping the M-theory circle. Notice 
that since CP 3 is a homogeneous space all positions for the IIA string are equivalent. 
The two computations (large N matrix model and area) of course agreen 

This Wilson loop is 1/2 BPS in a general M = 2 supersymmetric gauge theory on 
S 3 , as we review in section [2j and can be computed using the large N matrix model 
description. The supergravity dual computation will naturally involve an M2-brane 
wrapping the M-theory circle, leading to the same fundamental string configuration 
in Euclidean AdS4 (see Figure [I]). The only issue is which copy of the M-theory 
circle is relevant? When the internal space is Y? = S 7 /Z& all choices are equivalent 
by symmetry, but on a general Sasaki-Einstein manifold Yj, or a more general non- 
Einstein Y-j with flux, this is clearly not the case. Equivalently we may ask which IIA 
fundamental strings in AdS 4 x M 6 , that are pointlike in M e = Y 7 /U(1)m, preserve any 
supersymmetry. 

Summary of results 

Our starting point is to consider BPS M2-branes in general M = 2 supersymmetric 
AdS 4 x Y 7 solutions of eleven-dimensional supergravity. These backgrounds were stud- 
ied in detail in [24^ 25], where it was shown that provided the quantized M2-brane 
charge N of the background (measured by a certain flux integral) is non-zero, then 
there is always a canonical contact one- form rj defined on Y?. Concretely, 77 is con- 
structed as a bilinear in the Killing spinors on Yj, and it was shown in the latter 
reference that this contact structure entirely captures both the gravitational free en- 
ergy of the background, and also the scaling dimensions of BPS operators arising from 
supersymmetric M5-branes wrapped on five-manifolds £5 C Yj. 

In this paper we will show that the same contact form rj captures the Wilson loop 
VEV (W) of interest, computed semi-classically from the action of a BPS M2-brane. 



2 Similar Wilson loops have recently been considered in five-dimensional superconformal field the- 
ories on S 5 [IS], which may also be computed using localization techniques. The gravity duals are 
described by warped AdSg X S 4 /Z n solutions of massive IIA supergravity, and thus the geometry of 
the internal spaces here is fixed and in fact unique [44 . 



calibrated M-theory circles 




type IIA internal space 
M 6 = Y 7 /U(1) M 



supersymmetnc points 
in M 6 = fixed points of \ 



Euclidean AdS 4 







iU 




minimal surface 

2 2 




Hopf S'CS 3 



conformal boundary S 



Figure 1: A depiction of the total spacetime AdS4 x Y 7 , with a choice of M-theory circle 
U(1)m, together with the supersymmetric M2-branes of interest which are shown in 
red. These M2-branes are pointlike in the type IIA internal space M 6 = Y 7 /U(1)m, 
wrapping copies of the M-theory circle over these points, and are calibrated by the 
contact form rj. The supersymmetric points in M$ are precisely the points where 
the projection of the R-symmetry/Reeb vector field £ is zero (giving fixed points on 
M 6 ), and in general the calibrated circles over such points have different lengths. The 
remaining worldvolume of the M2-brane wraps a minimal supersymmetric surface S 2 
in Euclidean AdS4. The latter may be viewed as a hyperbolic 4-ball, with conformal 
boundary S 3 , and E 2 then has the topology of a 2-ball, with boundary a Hopf S 1 C S 3 . 



Concretely, we derive the general formula 

log (W) g rav it y = 



(27r) 2 / s i V , 

"" N 1/2 



V / 96VoT^ 
where we have defined the contact volume of Y-j as 



XD 



Vol,(y 7 ) = l^r ? A(dr ? ) 3 . 



(1.2) 



In particular, a supersymmetric M2-brane is calibrated with respect to r], which is why 



the integral of rj along the M-theory circle S X M appears in the formula ( 1.1 ). A contact 



form rj always has an associated unique Reeb vector field £, denned via the equations 
£j// = 1, £,-idi] = 0, and in [2U [25] it was shown that £ is also the R- symmetry 
Killing vector field, that is expected since an M = 2 superconformal theory in three 
dimensions has a u(l)# symmetry in the superconformal algebra. We will show that 
an M2-brane wrapping a copy of the M-theory circle S^ is supersymmetric precisely 
when the generating vector field (m of U(1)m is proportional to £. Geometrically, this 
means that the corresponding fundamental string at a point p G Me is supersymmetric 
precisely when p is a fixed point of £, considered as a vector field on M 6 (on Y 7 , on the 
other hand, £ is always nowhere zero). 

There is another way to describe which wrapped M2-branes are supersymmetric 
which involves the Hamiltonian function for the M-theory circle, defined as 

hu = v(Cm) • (1-3) 

This is a real function Hm : l^ — > M, invariant under (m, and we show that the 
supersymmetric M-theory circles S\ 4 C Y 7 lie precisely on the critical set dfiM = 0. 
The action of a supersymmetric M2-brane corresponding to a point p G M 6 may then 
also be written as 

S M2 = W h »® N V» t (L4) 

'96Vol„(Y 7 ) 



where p 6 Y? is any point that projects top 6 M 6 = y 7 /[/(l)M- Since (1.4) depends 



only on r\ we may compute this expression in examples using the same methods em- 



ployed in [231 [25], [45] [49]. For example, for toric solutions (1.4) may be computed 
entirely using toric geometry methods. In general there are multiple supersymmetric 
Sm circles, which can have different lengths with respect to rj and thus leading to dif- 



ferent actions (1.4). In the semi-classical computation one should sum over all such 



configurations, which in the large iV limit then implies that in (1.1) it is the longest 
Sm that gives the leading contribution to the Wilson loop. 

In the families of examples that we shall study, the dual field theory computation of 
the Wilson loop VEV reduces to a computation in a large N matrix model. As we shall 
review in section [2j in this matrix model the eigenvalues at large iV take the general 
form X 1 = xN 1 / 2 + iy 7 (x), where the index I runs over the number of factors of U(N) 
in the gauge group G = TJ 7 U(N), and are described by an eigenvalue density function 
p(x) which is supported on some interval [x min , x max ] C K. To leading order at large 
N it is straightforward to compute 

io g (iy) QFT = x max iv 1/2 , (i.5) 



which should be compared to the dual supergravity result (1.1) 



Remarkably, in all examples that we study we find that the interval [x m i n , x max ] in the 
matrix model coincides, in a precise way, with the image of the Hamiltonian function 
HmO^t)- Since Y-j is compact and connected, the latter image is also necessarily a 
closed interval, and more precisely we find tiM(Y 7 ) = [c min , c max ], where the field theory 
quantity x is proportional to the geometrical quantity c: 

x = , J =c. 1.6 

^eVol^Fr) 



The Hamiltonian Km is a Morse-Bott function on the symplectic cone over I7, and 
on general grounds we know that the image interval [c min , c max ] is divided into P 
subintervals c min = c\ < c 2 < ■ ■ ■ < c P+ i = c max , where the critical set maps as 
Km ({d/iM = 0}) = {ci I i = 1, .. .,P + 1}. For all c G (q,q + i) the level surfaces 
h^j(c) C Y-j are diffeomorphic to a fixed six- manifold, with the topology changing pre- 
cisely as one passes a critical point q. Even more remarkable is that we find that the 



corresponding points Xi, related to q via (1.6), are precisely the points where p'(x) has 



a jump discontinuity in the matrix model. These points are then also related to the 
fixed points of the Reeb vector £ on M 6 . 

The outline of the rest of this article is as follows. In section[2]we review the definition 
of the BPS Wilson loop in TV = 2 Chern-Simons-matter theories, and how it may be 
computed in the large N matrix model. Section [3] analyses supersymmetric M2-branes 
in a general class of AdS4 x Y-j backgrounds in M-theory, and we derive the general 



formula for the action (1.4), leading to the holographic Wilson loop result (1.1). In 
section [4] we compute the Wilson loop, on both sides of the correspondence, in a variety 
of examples, including for several infinite families of Sasaki-Einstein Y?, and for models 
with non-Einstein Y 7 . Section [5] concludes with a brief discussion. 

2 Wilson loops in J\f = 2 gauge theories on S* 3 

The dual superconformal field theories of interest have UV descriptions as M = 2 
Chern-Simons gauge theories with matter on S* 3 . We begin in this section by defining 
the BPS Wilson loop in such a theory, summarize how it localizes in the matrix model, 
and explain how it can be efficiently calculated. This section is mainly a review of 
material already in the literature. 

2.1 The Wilson loop 

In M = 2 supersymmetric gauge theories the gauge field A^ is part of a vector multiplet 
that also contains two real scalars a and D, that are auxiliary fields, and a two- 
component spinor A, all of which are in the adjoint representation of the gauge group 



dim 1Z 



(2.1) 



G. The BPS Wilson loop in a representation 1Z of G is given by 

Pexp ( (t ds(L4^ + a\x\) ) 

where x^(s) parametrizes the worldline 7 C 5* 3 of the Wilson line and the path ordering 
operator has been denoted by V. For a Chern-Simons theory the gauge multiplet has 
a kinetic term described by the supersymmetric Chern-Simons action 



^Chern-Simons = — / d 3 xy/detg Tr I A A dA + - A A A A A - A f A + 2Da J 



(2.2) 



where here g is the round metric on S 3 , and fc denotes the Chern-Simons coupling. 
When G is a product of unitary groups, G = YIjU(Ni), one can in general take 
different fcj G Z for each factor. In this case we will denote k = gcd{kj} [3]. 

There are four Killing spinors on S 3 , two satisfying each choice of sign in the equation 
V^e = ±|r^£, where the gamma matrices r M in an orthonormal frame generate the 
Clifford algebra Cliff (3, 0), and may thus be taken to be the Pauli matrices. A natural 
orthonormal frame {e m } m= i i 2,3 on S 3 is provided by the left (or right) invariant one- 
forms under the isometry group SU(2) Mt x SU(2) Tight . The four Killing spinors on S 3 
transform in the (2, 1), (1, 2) representations of this group. 

The full supersymmetry transformations for a vector multiplet and matter multiplet 
may be found in [261 E3 EE] ■ For our purposes we need note only that localization of 
the path integral, discussed in the next section, requires one to choose a Killing spinor 
e, which without loss of generality we assume solves V M e = |t m £- This choice of Killing 
spinor then has the two associated supersymmetry transformations 



SAfj. = -i^ r v e 



5a = -^e . (2.3) 



If one varies the Wilson loop (2.1 ) under the latter supersymmetry transformation one 
obtains 

6W oc \\\t^ - |i|)e . (2.4) 

The Wilson loop is then invariant under supersymmetry provided 

[r^ - \x\)e = . (2.5) 

Choosing s to parametrize arclength, so that \x\ = 1 along the loop, we see that t^x^ 
must be constant. In the left-invariant orthonormal frame e m one may then align x^ 



along one direction, say e 3 . The integral curve of this vector field is a Hopf S 1 C S 3 
(or equivalently a great circle). The supersymmetry condition then becomes 

(r 3 -l)e = 0. (2.6) 

This projection condition then fixes one of the two possible choices of e satisfying 



V^e = k r n £ t implying that the Wilson loop (2.1 ) is indeed a 1/2 BPS operator provided 



one takes 7 to be a Hopf circle. We will see later on that the condition (2.6), plus the 
fact that the supersymmetry generators are Killing spinors, also arises as the condition 
for supersymmetry of a probe M2-brane. 

2.2 Localization in the matrix model 



The VEV of the BPS Wilson loop (2.1) is, by definition, obtained by inserting W into 
the path integral for the theory on S 3 . The computation of this is greatly simplified 
by the fact that this path integral localizes onto supersymmetric configurations of 
fields. We summarize the main steps and results in this section, following in particular 
[271 EQl [31] , and refer the reader to the original papers for further details. 



The central idea is that the path integral, with W inserted, is invariant under the 



supersymmetry variation 5 corresponding to the Killing spinor e satisfying (2.6). We 



have written two of the supersymmetry variations in (2.3), and the variations of other 
fields (including fields in the chiral matter multiplets) may be found in the above ref- 
erences. Crucially, 5 2 = is nilpotent. There is then a form of fixed point theorem that 
implies that the only net contributions to this path integral come from field configu- 
rations that are invariant under S [SO]. Formally, one can argue this by introducing a 
collective Grassmann coordinate $ along the direction defined by 5 in field space, and 
then appeal to the fact that the Grassmann integral J di? = 0. This then breaks down 
precisely at fixed points of 5, where the coordinate d is not defined. 

Alternatively, and more practically for computation, one may add a conveniently 
chosen 5-exact positive definite term to the action, which a standard argument shows 
does not affect the expectation value of any supersymmetric (^-invariant) operator. For 
the vector multiplet one can add the term tTr[(5A)^5A] to the action (a similar term 
exists for a matter multiplet), without affecting the path integral. Sending t — > 00 one 
notes that, due to the form of this term added to the Lagrangian, only configurations 
with 5X = contribute to the path integral in a saddle point approximation. This 
saddle point then gives the same value as if the path integral had been calculated with 
t = 0, which is the quantity we are interested in. The saddle point approximation 
requires one to compute a one-loop determinant around the (^-invariant field configu- 
rations, which in the terminology of fixed point theorems is the contribution from the 
normal bundle to the fixed point set in field space. 

8 



For the M = 2 supersymmetric Chern- Simons-matter theories of interest, one finds 
that the 5-invariant configurations on S 3 are particularly simple: 



A,. 



, and D 



-a 



constant 



(2.7) 



with all fields in the matter multiplet set identically to zero. Here we may diago- 
nalize a by a gauge transformation. For a U(N) gauge group we may thus write 
a = diag(|r,... ^r), thus parametrizing 2na by its eigenvalues A,. The theories of 
interest will have a product gauge group of the form G = Yli=i U(N), and for t = oo 
the partition function then takes the saddle point form 




cxp 



where the one-loop determinant is given by 



k 



N 



i=i i=i 



one — loop 



(2.8) 



one — loop 



] | ] J 2 sinh 

i=i i^j 



M 



H 



]| detTz a exp[£(l-A a + ia)} . (2.9) 



matter a 



Here the first exponential term in (2.8) is simply the classical Chern-Simons action 



in (2.2), evaluated on the localized constant field configuration (2.7). The one-loop 



determinant factorizes, and the first term in (2.9) is the one-loop determinant for the 



vector multiplet. Since we have used gauge-invariance in (2.8) to restrict the integral 



to the Cartan subalgebra, we also have a Vandermonde determinant which has been 
cancelled against a term that appears in the one-loop determinant. The second term in 



(2.9) involves a product over chiral matter multiplets, labelled by a. We have taken the 



a multiplet to be in representation 1Z ai and with R-charge A a . The deteminant in 
the representation lZ a is understood to be a product over weights g in the weight-space 



decomposition of this representation, and a is then understood to mean g(a) in (2.9). 
Finally, 



£(z) 



■log(l 



Jl-KXZ 



K 



TTZ 2 + 



7T 



Li 2 (e 27ri2 ) 



17T 

12 



(2.10) 



In this set-up, the VEV of the BPS Wilson loop (2.1) reduces to 

g n 



(W) 



Z(N\)9 dimTZ 




e ^ s I=1 ^UH) 2 Tr n (e 2 ™)e 



-one — loop 



(2.11) 
Notice the integrand is the same as that for the partition function (2.8), with an 



additional insertion of Tr^e 27 ™") arising from the Wilson loop operator. The factor 



of (iV!) 5 , as in (2.8), arises from dividing by residual Weyl transformations, which for 



U(N) introduces a factor of 1/iV!. Note also that we have normalized the VEV relative 
to the partition function Z, so that (1) = 1, as is usual in quantum field theory. 
Localization has reduced the partition function Z and the Wilson loop VEV to finite- 



dimensional integrals (2.8), (2.11) over the eigenvalues Af of a, but in practice these 



are difficult to evaluate explicitly due to the complicated one-loop effective potential 



(2.9). For comparison to the dual supergravity results we must take the N — > oo limit, 
where the number of eigenvalues, and hence integrals, tends to infinity. One can then 
attempt to compute this limit using a saddle point approximation of the integral (this 
is then our second application of the saddle point method). With the exception of the 
M = 6 supersymmetric ABJM theory, where this matrix model is well-understood [51] , 
for general M = 2 theories the large N limit of the matrix integrals is not understood 
rigorously. However, in J30J a simple ansatz for the large N limit of the saddle point 
eigenvalue distribution was introduced. This ansatz is based on a partial analytic 
analysis of the matrix model, and also on a numerical approach to computing the 
saddle point. One seeks saddle points with eigenvalues of the form 

\\ = XiNP + iy!, (2.12) 

with Xi and y{ real and assumed to be 0(1) in a large N expansion, and j3 > 0. In the 
large N limit the real part is assumed to become dense. Ordering the eigenvalues so 
that the Xi are strictly increasing, the real part becomes a continuous variable x, with 
density p(x), while y{ becomes a continuous function of x, y\x). 



Substituting this ansatz into the partition function expression (2.8), the sums over 



eigenvalues become Riemann integrals over x, and one finds that the double sums 



appearing in the one-loop expression (2.9 ) effectively have a delta function contribution 
which reduces them to single integrals over x. (This is often described by saying that 
the long range forces in the matrix model cancel.) Writing Z = e~ F one then obtains a 
functional F[p(x),y I (x)], with x supported on some interval [x m i n ,x max ], and to apply 
the saddle point method one then extremizes F with respect to p(x), y T (x), subject to 
the constraint that p(x) is a density 

iax 

p(x)dx = 1 . (2.13) 



The existence of such a saddle point fixes the exponent (3 = | in (2.12). One then 
finally also extremizes over the choice of interval, by varying with respect to x m [ n , 
^max, to obtain the saddle point eigenvalue distribution p(x), y\x). 



We shall be interested in evaluating the Wilson loop VEV (2.11 ) in the fundamental 
representation, so that the Wilson loop is proportional to X]?=i S«=i eAi • I n the large 
N limit, described by the saddle point density p(x) and imaginary parts y\x) of the 

10 



eigenvalues, the VEV reduces simply to 

(^) QF t = NJ2 e xN/2+iyI Wp{x)dx. (2.14) 

T 1 *f 3Jmin 



1=1 



Because of the form of F[p(x),y I (x)} for M = 2 Chern-Simons-matter theories, the 
saddle point eigenvalue density p(x) is always a continuous, piecewise linear function 
on (x min , x max ). A simple computation then shows that, to leading order in the large 



N limit, the matrix model VEV (2.14) reduces to 



log(W% FT = x^N 1 ' 2 . (2.15) 



This is our final formula for the large N limit of the Wilson loop VEV. We see that it 
computes the maximum value of the (real part of the) saddle point eigenvalues. 

In our summary above we have suppressed the dependence on the R-charges A a of 



the matter multiplets, labelled by a, appearing in (2.9). If these are left arbitrary, one 
obtains a free energy F that is a function of A a , and according to [27] the supercon- 
formal R-symmetry of an M = 2 superconformal field theory further extremizes F as 
a function of A a (in fact maximizing F [52] )• For theories with M-theory duals of the 
form AdS 4 x Y 7 one finds the expected supergravity result 



2jV " -N 3 ' 2 , (2.16) 



27Vo\ v (Y 7 
but as a function of R-charges A a [31J, where on the right hand side it is in general 



the contact volume (1.2) of Y 7 that appears, as a function of the Reeb vector field £. 
This has by now been demonstrated in many classes of examples in the literature [25] , 

(HDH37J. 



3 BPS M2-branes 

In this section we analyse the supersymmetric probe M2-branes that are relevant for 



computing the holographic dual of the Wilson loop VEV (2.15). We first recast the 
condition of supersymmetry into a geometric condition, then derive the formula (1.4) 
for the action of the M2-brane, and finally describe how this may be computed in 
practice using different geometric methods. 

3.1 Supergravity backgrounds 

We will study the general class of M = 2 supersymmetric AdS 4 x Y? backgrounds of M- 
theory described in [241125] . We begin by recalling some relevant results and formulae. 

11 



The eleven-dimensional metric and four-form G 4 take the form 

9u = e2A ( j0AdS 4 + 9y t 

G 4 = -V0I4 + F4, (3.1) 

lo 

where the metric on AdS 4 here has unit AdS radius, with volume form V0I4. The warp 
factor A is taken to be a function on Y 7 , m is a constant, and _F 4 is a four-form on 
Y 7 . This is the most general ansatz compatible with the symmetries of AdS 4 . The 
eleven-dimensional Majorana spinor takes the form 

e = e ' 2 ip+ £g> x+ + e ip- <B> X- + charge conjugate , (3.2) 

where x± are complex spinors on Yj, ip± are the usual Killing spinors on AdS 4 (the 
± signs are related to the charge under the R-symmetry, discussed below), and the 
factors of e A / 2 have been introduced for convenience. 

In general the spinors x± solve quite a complicated system of coupled first order 
equations on Y 7 , that may be found in (2U |25]. These equations are then necessary 
and sufficient for supersymmetry of the AdS 4 x Y 7 background. For our purposes we 
need note only a few key formulae. We first define the real one-forms 

6 
f = ix+7(i)X- > V = e3A X+7(i)X+ , (3.3) 

where in general we denote 7(„) = ^ r y mi ... mn dy mi A • • • A dy 171 ", with y 1 , . . . ,y 7 local 
coordinates on Y 7 , and the superscript c on the spinors denotes charge conjugation. By 
an abuse of notation, we'll more generally regard £ as the dual vector field defined by 
the metric gy 7 - We then note that the differential equations for x± imply the equations 

777 

X+X+ = X-X- = 1 , "6" e ~ 3A = _Im [x+X-] » Re [x+X-] = ° , 
Re [x+7(i)X-] = ° > X+l(i)X+ = ~X-1{\)X- , 

dv = -^ 3A Re [xllmX-] ■ (3.4) 

These equations may all be found in reference [25] . 

The one-form 7] is a contact form on Y 7 , meaning that the top form r] A (drj) 3 is 
nowhere zero. Indeed, one finds 1251 that 



r] A(dr ] ) 3 = ^e 9A vol 7 , (3.5) 

where V0I7 is the Riemannian volume form defined by gy 7 - It is a general fact that a 
contact form 77 has associated to it a unique Reeb vector field, defined by the relations 

^r? = 1 , ^d V = , (3.6) 
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and remarkably one finds that £ and i] defined by (3.3) indeed satisfy these equations. 
Moreover, £ is a Killing vector field under which x± carry charges ±2, and as such is 
the expected R-symmetry vector field. 

Dirac quantization in this background implies that 

should be an integer, where £ p denotes the eleven-dimensional Planck length and G<± = 
dC*3. This may be identified with the M2-brane charge of the background, and a 
computation pH [25] gives 

N = ww^h^ f < (3 ' 8) 



relating the quantized M2-brane charge to the contact volume (1.2) of Y 7 and m. Since 
this is proportional to m 2 , in fact the contact form in (3.3) may be defined only when 



this charge is non-zero, so that m ^ 0. We assume this henceforth. 

The above supergravity solution of M-theory is valid only in the large N limit, even 
for solutions with non-trivial warp factor A and internal four-form flux F 4 . To see 
this [25] , note that the scaling symmetry of eleven-dimensional supergravity in which 
the metric g\\ and four-form G4 have weights two and three, respectively, leads to a 
symmetry in which one shifts A — y A + k and simultaneously scales m — y e 3K m, F 4 — y 
e 3K F, where k is any real constant. We may then take the metric gy 7 on Y 7 to be of order 



0(1) in N, and conclude from the quantization condition (3.8), which has weight 6 on 



the right hand side, and the expression for me _3A in (3.4) that e A = 0(N 1 ' 6 ). It follows 
that the AdS 4 radius, while dependent on Y 7 in general, is -RAds 4 = e A = 0(iV 1//6 ), and 
that the supergravity approximation we have been using is valid only in the iV — y 00 
limit. 

3.2 Choice of M-theory circle 

In addition to this background we must also pick a choice of M-theory circle. Geo- 
metrically, this means we also choose a U(l) = U(1)m action on Y7. At first sight it 
might seem to be contradictory that the supergravity computation we describe then 
manifestly depends on a choice of M-theory circle, while the dual superconformal field 
theory apparently does not. However, recall that the UV description of the gauge the- 
ory, whose Lagrangian we used to compute the localized path integral and Wilson loop 
in section 2.2, does in fact require a choice of M-theory circle U(1)m- We may have 
two or more such theories, arising from different choices of U(1)m and flowing to the 
same superconformal fixed point; but it does not follow that the Wilson loop operators 
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in these theories map to each other. One thus expects the Wilson loop VEV to depend 
on a choice of M-theory circle, in general. 

In terms of the supergravity solution described in the previous section, a choice of 
U(1)m implies the choice of a (to.oto.-U (1) r) Killing vector field (m on (Yj, <?y 7 ), whose 
flow generates the M-theory circle action. In particular Cm should preserve the Killing 
spinors x± on ^ii and hence also the contact one-form 77. The type IIA spacetime is 
then a warped product AdS4 x M 6 , where M 6 = Y 7 /U(1)m is the quotient space. 

Of course globally we must be careful when writing Mq = Yj/U(1)m- Although 
in principle one might choose any U(1)m action on Y?, in practice the gauge theories 
we study arise from "nice" actions of U(1)m- In particular, if the action is free then 
Mq inherits the structure of a smooth manifold from Y?, the simplest example being 
that of the ABJM theory with M 6 = CP 3 = S 7 /U(l) Uop {- If one embeds S 7 C C 4 
as a unit sphere in the obvious way, then recall that U(l)n op f may be taken to have 
weights (1,1,-1,-1) on the four complex coordinates (z±, z 2 , 23, 24) on C 4 . In this 
case the dual field theory is the M = 6 ABJM theory, which in M = 2 language is a 
U(N) x U(N) Chern-Simons gauge theory with two chiral matter fields Ai, A 2 in the 
bifundamental (N,N) representation of this gauge group, two chiral matter fields B\, 
B 2 in the conjugate (N, N) representation, and a quartic superpotential. 

Another important case is when U(1)m acts on Y-j with a codimension four fixed 
point set & C Y 7 , and is free on the complement of this fixed point set. In this case 
the action on the normal space M 4 to a fixed point is via (w x ,w 2 ) h-> (e ltp Wi,e lv W2), 
where locally (m = d v and (wi,w 2 ) are complex coordinates on IR 4 = C © C In this 
case the quotient normal space is IR 3 = M. 4 /U(1)m, with the fixed point set & at the 
origin becoming a D6-brane locus in the type IIA spacetime. With this understanding, 
the IIA spacetime is again a smooth AdS4 x M 6 , but with an AdS^filling D6-brane 
wrapping J^", now thought of as a submanifold & C Mq. Again, the simplest example 
is a quotient of S 7 , but now where U(1)m has weights (1, —1,0,0) on the coordinates 
(zi, z 2 ,Z3, Z4) on C 4 D S 7 . This fixes a copy of & = S 3 C S 7 at z x = z 2 = 0, which 
then becomes a D6-brane locus in the type IIA spacetime Mq = S 7 /U(1)m — S 6 . The 
dual field theory is then the low-energy gauge theory on iV D2-branes in flat spacetime, 
which in M = 2 language is a U(N) gauge theory with adjoint fields X ly X 2 ,X 3 and 
cubic superpotential (Af = 8 super- Yang-Mills) , but with additional fundamental fields 
arising from the low-energy string modes stretching between the D2-branes and D6- 
brane. This gives rise to additional fields q, q in the fundamental and anti-fundamental 
of U(N) respectively, and a corresponding additional superpotential term (see [T3]). 
This is often called the mirror to the ABJM theory, and indeed both theories have 
superconformal fixed points that are dual to AdS4 x S 7 . The gauge theories are of 
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course quite different, one being a U(N) x U(N) gauge theory, the other being a U(N) 
gauge theory]^] 

In the above cases the type IIA description is under control and typically well- 
understood, allowing one to determine an appropriate UV gauge theory. We shall see 
more complicated examples in section |4j 

3.3 BPS M2-brane probes 

The supersymmetric M2-brane which is conjectured to be holographically dual to the 
Wilson loop on S 3 must necessarily have as boundary a Hopf circle in S 3 . A convenient 
explicit form for the Euclidean AdS4 metric can be taken to be 

dq 2 
fi , Ads 4 = 1 2 + <? 2 d^3 , (3.9) 

with d^3 the round metric on the unit sphere S 3 , and q G [0, oo) a radial coordinate. 
The M2-branes of interest then wrap £ 2 x Sj^, where the surface £ 2 C AdS4 has 
boundary <9£ 2 = S^ opi C S 3 , and Sj^ C Yj is the M-theory circle. The submanifold 
£ 2 is then parametrized by the radial direction q in AdS,4, and a geodesic Hopf circle 
^Hopf m ^ 3 ' whilst S\ f C Yj is a priori arbitrary (imposing supersymmetry will later 
give restrictions on Sj/). The area of the surface £ 2 in AdS 4 is divergent, but can 
be regularized by subtracting the length of its boundary, i.e. the length of the S# opi 
geodesic in S 3 at q — > 00. Notice this is then a local boundary counterterm. Including 
also the warp factor one finds the regularized area to be 

Vol(E 2 ) = -^e 2A . (3.10) 

The action of the M2-brane then reads 

Vol(£, x SI,) 1 n f 3a 

*■ wn ~w?^hc vols: " (311) 



'M 



where V0I51 is the volume form on Sj^ induced from the metric gy 7 - 

As mentioned above, imposing that the M2-brane S 2 x Sj^ is supersymmetric gives 
restrictions on the possible circles S^. To see this, we need to split the Clifford al- 
gebra Cliff (11,0) generated by gamma matrices Ta satisfying {Ta,T b } = 26 ab into 
Cliff (4, 0) <8> Cliff (7, 0) via 

r Q = Pa®l, r a+3 = p 5 ®7a, (3.12) 



3 It happens that the Wilson loops turn out to be the same in these theories (essentially due to 
the high degree of symmetry), but the spectrum of BPS M2-branes/fundamental strings is certainly 



different. See section 4.1 
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where ot,/3 = 0, 1, 2, 3 and a,b = 1, . . . , 7 are orthonormal frame indices for Euclidean 
AdS 4 and Y 7 respectively, {pcnPp} = 2<5 a/ 3, {7 a ,7fc} = 2<5 afc and we have defined p 5 = 
PoPiP2P3- If we denote by X M the embedding coordinates of the worldvolume of the 
M2-brane into the target geometry, the amount of preserved supersymmetry is equal 



to the number of spinors e, as in (3.2), satisfying the projection condition 



Pe = 0, where P = ±(l-^e ijk d i X M d j X N d k X p r MN p\ , (3.13) 

with i, j, k indices on the worldvolume. We now choose an orthonormal frame in eleven- 



dimensions as (c.f. (3.1 



2 vTT? 2 



E° = ^e A , q , E m = ^e A qe m , E 3+a = e A e% , (3.14) 



where {e m } m= i i 2,3 is an orthonormal frame on S 3 and {ey .} a =i,...7 is an orthonormal 
frame on (Y 7 ,g 7 ), with ey 7 (or rather its dual vector field) aligned along the M-theory 



circle vector field Cm- Taking e 3 to be aligned along the Hopf circle, as in section 2A_ 
the projector P then takes the simple form 

P = ^(l-ip 5 P03®7i) , (3.15) 

and the constraints that follow on the spinors i/)±, x± on Euclidean AdS 4 and Y 7 , 
respectively, are 

(1 - ip 5 po3)ip± = 0, and (1 - 7i)x± = 0. (3.16) 

In order to determine how much supersymmetry is preserved by the brane in AdS4, 
we must count the number of Killing spinors ip± that satisfy the last projection equa- 
tion. We may decompose the four- dimensional gamma matrices into po — 1 ® T 3 and 
Pfj, = Tfj, <S> Ti, with the Pauli matrices 77, p = 1, 2, 3. These matrices act on spinors of 
the form ip = (V'i,'02) T , with ^ 1)2 2-component spinors. The Killing spinors on AdS 4 
may then be constructed from Killing spinors on the S 3 at fixed radial coordinate q. 
Explicitly, if e solves the Killing spinor equation 



on S 3 , then 



V^ = -V-, (3.17) 



"= '(^v^T?)-^! (3 ' 18) 



is a Killing spinor on Euclidean AdS 4 . Equation (3.17) has two solutions, one being 



chiral and one anti-chiral, i.e. t^e = ±e. One then easily shows that the first projection 



equation in (3.16) is satisfied if we restrict to chiral e in the last solution, which singles 
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out one of these two spinors on AdS^ Hence the M2-brane preserves half of the 
supersymmetry in ACIS4. Note that the same positive chirality condition also appeared 
in the supersymmetry condition derived in the field theory context, c.f. (2.6[)j^ 



The second projection equation in (3.16) tells us which circles S\j give rise to 



supersymmetry-preserving M2-branes. Following a standard argument one notices that 



X+ 



1 -71 



Xa 



Xa 



7i 



7i 



X-\ 



1 -7! 



Ah 



> , (3.19) 



using 7! = 7 J and 7^ = 1. This immediately gives vol 5 i > x+l(i)X+ (with a pull-back 
understood), with equality if and only if some supersymmetry is preserved by S 1 ^. The 
action ( |3.11 ) for a supersymmetric brane is then 

Voi(s 2 x s x M ) : 



7T 



>M2 



(27T) a fl 



(2tt) 2 ^ 2 J sl 



e 3A X+7(i)Xn 



(3.20) 



With the help of equations (3.3) and (3.8) the action of a supersymmetric M2-brane 



can be rewritten in terms of the contact form 77 as (taking a convention in which m < 0) 



•>M2 



2j Y7 r ] A(dr 1 y 



-.N 1 ' 2 . 



(3.21) 



3.4 M-theory Hamiltonian function 

In this subsection we further elucidate the geometry associated to these supersymmetric 
M2-branes. This geometric structure will both be of practical use, when we come to 



compute the M2-brane actions (3.21) in examples, and also, as we will see, is realized 



rather directly in the large iV dual matrix model. 

We begin by introducing the M-theory Hamiltonian function 



h 



M 



= v(C 



M, 



(m->V 



(3.22) 



where (m generates the M-theory circle action. This is a real function on Y 7 , and 
since (m is assumed to preserve the Killing spinors and metric on Y7, it follows that 
(ju preserves Hm and commutes with the Reeb vector field £. It follows that the 
contact length of an M-theory circle S]^ over a point p e M 6 = Y 7 /U(1)m is given by 
f s± r\ = 27r/i M (p), where p € Y 7 is any lift of the point p. This directly leads to the 



form of the M2-brane action (1.4). 



4 The other two Killing spinors on AdS4 are constructed from spinors on S 3 satisfying V^e = — ^r M e. 
We set the corresponding spinors to zero in section [21 as they are not used in the supersymmetric 
localization. Again, one chirality is broken by the M2-brane. 

5 Notice that the Wilson loop circle 7 C S 3 is calibrated by e 3 , one of the left-invariant one- forms 
under SU(2)\ c [ t , which is a contact form on S 3 . 
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One way to characterize the supersymmetric M-theory circles S M is to note that on 
TY 7 \ s i the vector Cm is necessarily proportional to the Reeb vector. Indeed, using 



(3.4) one can show that at these supersymmetric points 

(M^dri = 0. (3.23) 



To see this one takes the projection condition (3.16) with x—> applies x+la on the 



left, and then takes the real part of the resulting equation. Using Re [x+X- 



Re [x+7aX-] — and the relation between dr] and Re [x+7(2)X-] m (3-4) then leads 



to (3.23). That this then implies Cm oc £ follows from the fact that rj is a contact 
form: dt] is a symplectic form on ker?7, the rank 6 subbundle of the tangent bundle 
TY-j of Y-j defined as vectors having zero contraction with rj. Since this means that drj 
is non-degenerate on this rank 6 bundle, and since also TYj = kerrj © (£), where (£) 



is the real line bundle spanned by vectors proportional to £, (3.23) implies that the 
projection of Cm onto kerrj is zero, i.e. that Cm oc £. 



The condition (3.23) is then also the condition that we are at a critical point of the 



Hamiltonian Hm- To see this, we rewrite C^ M r] = using the Cartan formula, so that 



(3.23) is equivalent to 

d(Cj*J7?) = & dh M = . (3.24) 

Thus the supersymmetric M2-branes lie precisely on the critical set {dhM — 0}, and 



their action (1.4) is determined by Hm evaluated at the critical point! It is a general fact 
that any component of the moment map for a compact group action on a symplectic 
manifold is a Morse-Bott function. Here more precisely recall that the cone C(Y) = 
M>o x Y-j is symplectic, with symplectic form 

co = l -d(r 2 r,) , (3.25) 

where r > is a radial coordinate. In fact the cone being symplectic is equivalent to 
(Y 7 , rj) being contact. The M-theory circle action then gives a U(1)m action on this 
cone, with moment map 

fi = \r 2 Cu^. (3.26) 

Thus fj, is Morse-Bott, and the restriction of // to Y? at r = 1 is our Hamiltonian 
function /im/2. We thus know that the image /i m (Yt>) = [c min , c max ] is a closed interval, 
and this is further subdivided into P intervals via c min = c\ < c 2 < ■ ■ ■ < c P+1 = c max , 
where the q are images under Km of the critical set {dhM = 0}. On each open 
interval c G (cj,q + i) the level surfaces h M (c) are all diffeomorphic to the same fixed 
six- manifold, with the topology changing as one crosses a critical point q. 
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Finally, since at a supersymmetric S\f we have (m oc £, it follows that the corre- 
sponding point p G Mq = Y 7 /U{1)m is a fixed point under the induced Reeb vector 
action on M 6 = Y 7 /U(1)m- That is, over every fixed point p G M 6 of £, there exists 
a calibrated and supersymmetric M-theory circle S\j whose corresponding supersym- 



metric M2-brane action is given by (1.4). 



In the holographic computation of the Wilson loop VEV via the M2-brane action, 
one should sum e~ Sm2 ' p over all contributions. In some cases (typically with more 
symmetry) we shall find that the supersymmetric points p G Mq form submanifolds 
which are fixed by £, and this sum in fact becomes an integral over the different 
connected submanifolds. Notice that Hm is constant on each connected component 
of the fixed point set. In any case, in the large N limit only the longest circle S\,j 
survives, the others being exponentially suppressed relative to it in the sum/integral, 



hence proving formula (1.1) 



The calculation of the action of a supersymmetric M2-brane can be completely car- 
ried out once the Reeb vector field £ and the M-theory circle generator (m are known. 
Indeed, the contact volume Vo\. v (Y 7 ) is a function only of the Reeb vector [45j, and 
the length of a calibrated circle f ql r\ = l-Khuip) depends only on £, ( M and the 

M.p 

point p. Even though this could appear to be involved, the computation of these two 
quantities is relatively straightforward for appropriate classes of Y 7 . In particular, if we 
focus on toric manifolds, standard geometrical techniques can be exploited to straight- 
forwardly find all calibrated circles, i.e. the connected components of the critical set 
{d/ijvf = 0} C Y 7 , as well as the contact volume jl9]. This is the subject of the next 
subsection. 

3.5 Geometric methods of computation 

In this section we explain how to compute the various quantities we have been dis- 
cussing in appropriate classes of examples. We focus our discussion on toric geometries, 
which means that U(l) 4 acts on Y 7 , preserving the contact form rj. In this case there 
are some pretty geometric methods, first developed in J1B1IIH], that may be utilized to 
calculate the length of the calibrated M-theory circles, as well as the volumes of the 
internal spaces. We will thus focus on this class of solutions, although we note that 
the more general methods described in [49j may be used to attack non-toric cases. 

Let us begin with the symplectic cone (C(Y) = IR> x Y, u> = ld(r 2 r/)), but in general 
dimension In. Equivalently, (Y,i]) is contact with dimF = 2n—l. The toric condition 
means that U(l) n acts on the symplectic cone C(Y) preserving the symplectic form 
u, and we may parametrize the generating vector fields as <9<^, with <pi G [0, 27r) and 
i = 1, . . . ,n. This allows one to introduce symplectic coordinates (yt, (pi) in which the 
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symplectic form on C(Y) has the simple expression 

n 

^d^Ad^. (3.27) 






Moreover, when the toric cone is of Reeb type, meaning that £ is in the Lie algebra 
of U(l) n , the coordinates y, take values in a convex polyhedral cone C* C M. n |54j . 
If this cone has d facets, we have corresponding outward primitive normal vectors to 
these facets, v a e Z n , a = 1, . . . , d, with the facets corresponding (under the moment 
map) to the fixed point sets of U(l) C U(l) n with weights v a . In particular this set-up 
applies to toric Sasakian Y [48J, in which the symplectic cone C(Y) is also Kdhler. In 
this case, the topological condition that C(Y) is Calabi-Yau (more precisely, that the 
apex {r = 0} is a Gorenstein singularity) is equivalent to the existence of an SL(n,Z) 
transformation such that the normal vectors take the form v a = (l,w a ), for all a, with 
w a G Z n_1 . In this basis, the first component of the Reeb vector is necessarily £i = n 



In general the components of £ = Y^i=x^<h f° rm a vector £ = (£i,... , £ n ) that 
defines the characteristic hyperplane in M. n : {y e M n | £ • y = |}. This hyperplane 
interesects C* to form a finite polytope A^, and the contact volume of the base Y is 
related to the volume of this polytope bj 



Vol„(y) = 2n(2vr) n Vol(A 5 ) . (3.28) 

Moreover, each of the d facets J-" a , intersected with the characteristic hyperplane, are 
images under the moment map of (2n — 3)-dimensional subspaces E a of Y. The volumes 
of these submanifolds may be calculated once the volumes of the facets are known, for 

Vol„(E a ) = (2n-2)(27rr- 1 ^Vol(J- a ) . (3.29) 

In addition, the volume of the base manifold Y is simply given by 

VoL,(F) = i^ ]T ^Vol( j;) . (3.30) 

In [39] the idea is to study the space of Kahler cone metrics on C(Y), and thus 
Sasakian structures on Y. One then considers the Einstein-Hilbert action (with a fixed 
positive cosmological constant) restricted to this space of Sasakian metrics on Y, so 
that a Sasaki-Einstein metric on Y is a critical point. In fact the action is minimized 
and proportional to the volume of the base Vol(Y) when the metric on Y is Sasaki- 
Einstein. In this case there is unique Reeb vector of the form £ = (n,^2, ■ ■ ■ , £n) such 



In the Sasakian case the Riemannian volume and contact volumes coincide. 

20 



that the Einstein-Hilbert action, or equivalently Vol(Y), is minimized as a function 



of £. Thus, for any given toric diagram one calculates Vol(Y) with formula (3.30) as 
a function of the Reeb vector, and determines £ for the Sasaki-Einstein metric on Y 
by minimizing this function]^] Presumably these ideas extend to more general warped 



geometries, with non-zero internal flux F4 7^ in (3.1), following a similar construction 
in type I IB AdSs solutions [45] . 

In this paper we need only apply this method for n = 4. A way to compute Vol(J-' a ) 
as a function of the Reeb vector for n = 4 has been described in [8]. If the facet 
J- a is a tetrahedron, its vertex is at the origin in C* and its base is a triangle lying 
in the characteristic hyperplane. This is generated by three edges passing from the 
characteristic hyperplane to the origin, and bounded by four hyperplanes creating the 
polyhedron. In addition to v a , three other facets are then involved in the construction 
of the tetrahedron, and we denote their normal vectors as u ttj i, v a> 2, v a ^. The volume of 
the tetrahedron may be expressed as 

1 VUT , \ l (Va,Va,l,Va,2,Va,3) 2 ,„ „■, \ 

\Va\ 48 I (£, V a , V atl ,V at 2){€, V a , V a ,l, Uo,3)(?, V a , V a>2 , V a ,3) \ 

with (•,-,•,•) the determinant of a 4 x 4 matrix. If the facet J- a is not a tetrahedron, i. e. 
there are more than 3 edges that meet at a vertex in the toric diagram (c.f. below), 
the volume can be computed with the same formula by breaking up the facet into 
tetrahedrons. 

The toric diagram for a toric Calabi-Yau cone is by definition the convex hull of the 
lattice vectors w a in n — 1 =3 dimensions. To each vertex in this diagram corresponds a 
facet J- a . If the vertex is located at the intersection of three planes, or equivalently three 
edges of the toric diagram meet at the vertex, then it corresponds to a tetrahedron. If 
instead four edges meet at the vertex, the facet is a pyramid than can be split into two 
tetrahedrons, and so on. A given facet T a then corresponds to a vector v a = (l,w a ), 
with w a a vertex in the toric diagram; the other three vectors v a; i,v a> 2,v a> 3 are the 
outward-pointing primitive vectors corresponding in the toric diagram to the three 
edges that meet at the vertex v a . Let us also note that the base Y 7 of the cone is a 
smooth manifold only if each face of the toric diagram is a triangle, and there are no 
lattice points internal to any edge or face. These conditions are equivalent [56] to the 
cone being good, in the sense of [54] . 

It should now be clear that once a toric diagram is given for a toric Calabi-Yau 
cone C(Y), one can calculate the volume of the base Vol^(Y) as a function of the toric 
data and the Reeb vector that is parametrized by £ = (4, £2, £3, £4)- After minimizing 
the volume with respect to £2, £3, £4, one obtains the Reeb vector and Vol^(Y) as a 



7 That the Sasaki-Einstein metric indeed always exists was proven in 
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function of the toric data only. For more general warped solutions with flux the cone 
is not Ricci-flat Kahler, but in the examples we shall study later in section [4] the Reeb 
vector £ and toric contact structure are in fact known [25] . 

Next we turn to the M-theory Hamiltonian function Hm-, and the computation of 
the calibrated circles in Y 7 and their lengths. This involves, by definition, the choice 



of an M-theory circle Cm, as described in section 3.2 As we proved in the last section 



supersymmetric calibrated Sj^ exist where Cm is parallel to £. This is equivalent to 

Cm = v((m)Z = h M i , (3.32) 

as follows by taking the contraction of each side with rj. We can conclude that if we 
know the proportionality constant between Cm and £, the length of the corresponding 
calibrated M-theory circle, located over a fixed point p under £ in M 6 , is then simply 
27t/im(p) with p £ Y? any point projecting to p. In terms of the toric geometry above, 
notice that 

n 

h M = 2j2viCu, (3-33) 

where Cm = YH=i Cm^»- This may be regarded as a function on the polytope A^, that 
is the image of Y 7 under the moment map. 

The only remaining question is how to find where the two vectors Cm, £ are propor- 
tional to each other, or equivalently what are the critical points of Hm-, and also what 
is the value of Hm at those points. With the formalism at hand, this is straightforward 
to answer. Once a toric diagram and Cm are given, the Reeb vector (and the volume) 
can be found with the method described above. We may then find the solutions to the 
equation 

Cm = /3£ + ^a a w a , (3.34) 

with j3, a a real numbers, and I C {1, . . . , d} a subset of three facets which intersect. 
Geometrically, the intersection of three facets defines an edge of C*, which corresponds 
to a circle S 1 C Y 7 . This circle is a fixed point set of U(l) 3 C U(l) 4 defined by the 
three vectors v a , a e /, meaning that the generating U(l) vector fields corresponding 
to v a are zero over this circle, and hence Cm is parallel to £. Thus this 5* 1 is precisely a 
calibrated circle. The proportionality constant is then hu = v((m) = (3, and its length 
is 27t/im- Thus our problem boils down to linear algebra on the polyhedral cone. 



We make a few further geometrical observations. First, if (3.34) holds with /3 — 
then Cm actually fixes the S 1 , meaning that there must be D6-branes present. The 
M-theory circle then has zero length on such loci, formally leading to M2-branes with 



zero action; if Cm acts freely on Y 7 this cannot happen. Next we note that (3.34) cannot 
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hold with a a = for all a E I, since then ( M would be everywhere parallel to £, and 



this cannot happen since Cm is a non-R symmetry. However, it may happen that (3.34) 
holds with one or two (but not all three) of the coefficients a a = 0. Geometrically, 
this means that in this case (m is parallel to £ over the intersection of (respectively) 
the corresponding two or one facets with non-zero a a coefficients, leading to three- 
dimensional or five- dimensional subspaces of Yj which are fibred by calibrated S^ 
circles. These then descend to two-dimensional or four-dimensional fixed point sets of 
£ on M 6 = Y 7 /U(1)m, respectively. We shall see examples of this in the next section. 
Finally, if the toric diagram contains faces which have more than three sides, then 



(3.34) may hold for / being the corresponding set of 4 or more vectors v a . In this case 
the manifold has a locus of (worse than orbifold) singularities along the corresponding 
S 1 in Y?, and our theory above doesn't directly apply to these singular circles (their 
tangent spaces are not even a quotient of IR 7 ). 

Even though the above theoretical background may appear cumbersome, it is effec- 
tively not difficult to find the volume of Y?, its Reeb vector £ and all the calibrated 



circles and their lengths. Thanks to equation (1.4), the action for each corresponding 
M2-brane follows straightforwardly, and can be compared to the data extracted from 
the matrix model of the dual field theory. We examine these computations, in a variety 
of examples, in the next section. 

4 Examples 

In this section we illustrate the duality between geometries and matrix models in a 
wide variety of examples. In particular we will compute the image of the M-theory 
Hamiltonian /im(^) = [cmin, c max ], and show that it coincides with the support of the 



matrix model eigenvalues [x min ,x max ] via (1.6). The critical points of h M will be shown 



to map to the points x = Xi where p'{x) has a jump discontinuity, with the matching 
of Wilson loops being a corollary of this result for x = x E 



max- 



4.1 Duals to the round S 7 

We begin by studying two superconformal duals to AdS 4 x S 7 , where S 7 is equipped 
with its standard round Einstein metric. These differ in the choice of M-theory circle 



U(1)m acting on S 7 , as we discussed briefly in section 3.2. In this case the geometry is 
particularly simple, allowing us to illustrate the geometric structures we have described 
very explicitly. 
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4.1.1 ABJM theory 

The ABJM theory [TJ is an M = 6 superconformal U(N)k x U(N)^ k Chern-Simons- 
matter theory. In M = 2 language, there are two chiral matter fields A\, A 2 in the 
bifundamental (N,N) representation of this gauge group, two chiral matter fields B±, 
B 2 in the conjugate (N, N) representation, and a quartic superpotential. Here the 
subscript k G Z in U(N) k denotes the Chern-Simons level for the particular copy of 



U{N), as in Q. This theory is dual to AdS 4 x S 7 /Z k with N units of flux g§ , 
where Z fc C C/(l) Hop f = ^(1)m- 

We may realize S* 7 as the unit sphere S C M 8 — C 4 and take C/(l)nopf to have 
weights (1,1,-1,-1) on the four complex coordinates (zi, z 2 , z 3 , z 4 ) on C 4 . In this 
description the U(1)r symmetry of the Af = 2 subalgebra of the TV = 6 manifest 
superconformal symmetry of the theory has weights (1, 1, 1, 1) on C 4 , which gives an 
different Hopf action on C 4 . In [lj the variables C a = (A-]_, A 2 , B\, B%) were also used. 
In this choice of complex structure on IR 8 = C 4 the U(1)m and U(1)r weights above 
are interchanged; in these variables the SU(4)r symmetry of the theory, which acts 
isometrically on CP = {S 7 C C 4 }/£/(1)m, is manifest. However, to be uniform with 
the other examples we shall study, we shall fix the first complex structure on IR 8 = C 4 
above. 

In these coordinates S 7 = {(z±, z 2 , z 3 , z&) G C 4 | |-2i| 2 + \z 2 \ 2 + | ^^3 1 2 — I— \z^\ 2 = 1}, while 
the M-theory Hamiltonian function on 5* 7 /Zfc is 

h M = \{\zi\ 2 + \z^-\z^-\ Zi \ 2 ) . (4.1) 



In the toric geometry language of section |3.5[ we have the symplectic coordinates 
Ui — \\ z i\ 2 - The level sets h~^(c) are diffeomorphic to S* 3 x S 3 /Z k for c G (— |, |). 
Perhaps the easiest way to explain this is to note that dividing the levels sets also by 
U(1)m gives the Kahler quotient description of T 1,1 ] 8 ^ The level sets are then a circle 
bundle over T 1 ' 1 ^ S 2 x S 3 , with first Chern class keZ = H 2 (S 2 x S 3 , Z), which means 
they are diffeomorphic to S 3 x S 3 /Z k . Notice that these level sets are also described 
by 

\zi\ 2 + \z 2 \ 2 = i(l + cfc), \z 3 \ 2 + \zi\ 2 = ^(1 - ck) . (4.2) 

When c — > ±| the S 3 x S 3 /Z k level sets thus collapse to two copies of S 3 /Z k at 
{2:3 = z 4 = 0} and {z\ = z 2 = 0}, respectively. Thus the image h M (S 7 ) = [— |, 1] , 
with the endpoints c max = — c min = \ being the only two critical points of the Morse- 
Bott function Km- 



8 Of course this is directly related to the construction of the ABJM theory itself, as the M-theory 
lift of the theory on N D2-branes at the conifold singularity C(T 1,:L ), with k units of RR two-form 
flux through the vanishing S 2 . 
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The contact form in these coordinates is 

4 



E 



2 _ \^ I |2 

r = > k 



(4.3) 



Being Einstein, the contact volume of S' 7 /Z/ c is equal to the Riemannian volume, with 



Vo\(S 7 /Z k 



71 

3k 



(4.4) 



Our general formula (1.6) thus implies that the matrix model variable x should be 



related to the geometric quantity c above via 

(2vr) 3 



x 



^96Vo\(S 7 /Z k 



c = Txv2kc 



(4.5) 



The large N saddle point eigenvalue distribution for the ABJM theory was given in 
. The eigenvalues for the two gauge groups are related by 



X\x) = X 2 (x) = xN 1/2 + iy(x) 



where 



p{x) 



Vk 



y{x) 



Vk 



2ttV2 ' 2^ 

and the eigenvalues are supported on [x min , x max ], where x ma 

P(x) k 



x 



"Xr 



7T 



(4.6) 

(4.7) 
y/2/k. This 



x.. 



-^-n 



Figure 2: Eigenvalue density p(x) for the ABJM theory. 



of course agrees with the geometric formula (4.5), and since the density p(x) is constant 
on (x m ; n , x max ) (see Figure |2| its derivative is in particular continuous on this region. 



It is then automatic that the gravity formula (1.1 ) agrees with the field theory formula 



(1.5) for the Wilson loop, giving in both cases 



log (W) = TT\/^N 1/2 



(4.8) 
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4.1.2 Mirror theory 



As discussed in section 3.2, the mirror to the ABJM theory (with k — 1) arises by 
choosing a different M-theory circle action on S 7 . The field theory [13] is M = 8 U(N) 
super- Yang-Mills theory coupled to two additional fields q, q in the fundamental and 
anti-fundamental representation of U(N), respectively. The superpotential is 

W = Ti(qX 1 q + X 3 [X 1 ,X 2 }) , (4.9) 

where Xi, X 2 , X 3 are the adjoint chiral fields of the M = 8 theory, in M = 2 language. 

In this case the M-theory circle U(1)m has weights (1, —1,0,0) on S 7 C C 4 , which 
has a codimension four fixed point set & = S 3 = {zi = z 2 = 0} C S 7 . It follows that 
the type IIA internal space is Mq = S 6 , with a space-filling D6-brane wrapping a copy 
of S 3 C S 6 . The field theory described in the previous paragraph is then the theory 
on N D2-branes in flat space (A/" = 8 super- Yang-Mills) , but coupled to additional 
massless fields q, q arising from the lowest excitations of strings stretching between the 
D2-branes and D6-brane. 

Although the background geometry is the same as in the previous subsection, the 
M-theory Hamiltonian is nov^ 



h M = \ Zl \ 2 -\z 2 \ 2 . (4.10) 

The level surfaces h^(c) are described by 

2\ Zl \ 2 + \z 3 \ 2 + \z±\ 2 = 1 + c, 2\z 2 \ 2 + \z 3 \ 2 + \z±\ 2 = 1-c, (4.11) 
so that c £ [—1, 1]. However, the critical point set of Hm is quite different to that for 



the ABJM model (4.1). The endpoints c = +1, c = —1 are now the copies of S 1 C S 7 
at {z 2 = z 3 = Zi = 0} and {z\ = z 3 = z<± = 0}, respectively. (Compare to the ABJM 
model, where for k = 1 also c £ [—1, 1], but with the endpoints being images of copies 
of S 3 , rather than S 1 .) Moreover, there is an additional critical point at c = 0, which 
then intersects the D6-brane locus/fixed point set at {zi = z 2 = 0}. Indeed, on S 7 we 
have 

dh M = (zidzi + Zi&zx) - (z 2 dz 2 + z 2 dz 2 ) , 
4 
= ^2( Zi dzi + Zidzi) <=> = dr . (4.12) 

9 We could similarly choose to quotient by Z^ C U(1)m, which would lead to k D6-branes in the 
type IIA description. However, here we restricted to k = 1 in order to compare to the k — 1 ABJM 
theory, which is also dual to AdS4 x S 7 (the point being that the Z& quotients in each case are 



different). In fact the general k case is a = k, b = of section 4.3 
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Thus in addition to the endpoints {z 2 = z 3 = z 4 = 0} and {zi = z 3 = z± = 0}, we also 
have dh M = at {z\ = z 2 = 0} = S 3 , which is the fixed point set of U{1)m where 
hu = 0. Thus we have the three critical points c\ = c m ; n = — 1, c 2 — 0, c% — c max = 1. 
The topology of the level sets h~j^(c) is the same for c G (—1,0) and c G (0, 1), but 
with different circles collapsing on each side. For c G (0, 1) we may "solve" hu = c as 



z i 



\z 2 \ + c > 0, and note that consequently Z\ ^ on this locus. From (4.11) it 



follows that h A j(c) = S\ x 5* 5 , where S\ is parametrized by the phase of z\ 
On the other hand, for c G (—1, 0) instead we solve hu — c as \z 2 \ Ji 



\Zi\e v 



\zi\ 2 — c > 0, so 



that h M (c) 



5*2 x S 5 , where S^ is parametrized by the phase of z 2 = \z 2 \e 



id>2 



The general formula (1.6) implies that the matrix model variable x should be related 



to the geometric quantity c again via 



x 



(2*0 
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\/2« 



^96Vo\(S 7 ) 
which is the same formula as for the ABJM model with k 



(4.13) 
1. The large N saddle 



point eigenvalue distribution is in fact a special case of the models in section 4.3, with 
a = 1, b = in the notation of that section, and appears in [33]. In this case there is 
only a single gauge group, and one finds the eigenvalue density 



2n 2 \-*^ ^1111117 



where Xr, 



~ dj\~\ 



7T\/2, thus agreeing with (4.13). Moreover, the derivative of p 
p(x) A 




Figure 3: Eigenvalue density as a function of x. There are three points where p'(x) is 
discontinuous, corresponding to critical points of hyi- 

is discontinuous at the endpoints and at the point x = 0, which arises as the image of 



the D6-brane locus under hu- The Wilson loop is again given by (4.8), with k = 1. 



4.2 Dual to Q w /Z k 

Our next example is that of the homogeneous (and toric) Sasaki-Einstein manifold 
Q 1 ' 1 ' 1 /Zfe. The manifold Q 1 ' 1 ' 1 is the total space of an S* 1 fibration over the product 
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of three copies of S 2 , i.e. S 1 <-$■ Q 1,1,1 — ^ S 2 x S 2 x S 2 , which describes its structure 
as a regular Sasaki-Einstein manifold. Even though this manifold is toric, and the 
geometrical techniques described in section |3.5| can be applied, we will instead take 
advantage of the fact that the metric is known explicitly on this space. 
The Sasaki-Einstein metric on Q 1 ' 1,1 can be written as 

(3 \ 2 3 

dV + ^cosfcd^] +^(d0? + siii 2 i d^), (4.15) 

8=1 / 1=1 

where the coordinates 9i G [0, n] and ifi G [0, 2n) are the usual S 2 coordinates, and the 
coordinate ip G [0, An) parametrizes the S l fibre. The contact form is simply 

^cos^d^j , (4.16) 

and for the field theory model below the M-theory circle is generated by (m = z(9 vi + 
dy 2 ). The M-theory Hamiltonian follows straightforwardly and reads 

h M = v((m) = ^(cos#i + cos# 2 ) . (4.17) 

The length of a supersymmetric M-theory circle is always given by 27r/iju(p), where 
p G Y 7 covers a fixed point p of £, with p G M 6 = Y 7 /U(1)m- However, when the 
Sasaki-Einstein manifold is regular, as in the case at hand, we may also describe the 
supersymmetric M-theory circles in terms of the base Kahler-Einstein manifold Bq = 
Yj/U(l)n, where U(1)r is generated by the Reeb vector £. In this point of view, the 
supersymmetric M-theory circles cover fixed points of (m on B 6 , which in the case at 
hand is Bq = S 2 x S 2 x S 2 because £ = 4d^. These points are located at {(#i, 62) \ Q\ G 
{0,7r}, 9 2 G {0,7r}}. Thus one obtains three critical values c\ = c m j n = —-hi, C2 = 0, 
C3 = c max = =-r. Notice these are S 2 loci of critical points, parametrized by (6*3, <^ 3 ). 

Being Einstein, the contact volume of Q 1 ' 1,1 /^ is equal to the Riemannian volume, 
with 

VaKQW/Z*) = ^ , (4.18) 

and as usual the Z fc quotient is along U(1)m generated by Cm- The general formula 



(1.6) tells us that the matrix model variable x max = — a^min predicted from the gravity 



calculation is 

x = W c = 2^ (419) 

max v/96Vol(QM,i/ Zfc ) Cmax - V3* ' [ ' 

A dual field theory to Q 1 ' 1 ' 1 /^ has been proposed in [HI [T3]. This theory is closely 
related to the ABJM theory. In addition to the bifundamental fields A iy B iy a pair of 
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field in the (anti-) fundamental representation is added to each gauge group node, and 
one adds a cubic term to the superpotential 



>V cub ic = Tr (q 1 A 1 q 1 + q 2 A 2 q 2 ) 



(4.20) 



The corresponding matrix model has been worked out in [32] , where it was found that 
the density of the real part of the eigenvalues is 

k 



p(x) 



4 f 2"( 2Xmax 



•Xj J L\JL 1L1 J^QJJ^ "\ **J "\ **J Yi 



(4.21) 



with x r 



2tt 



3k 



, thus agreeing with (4.19). Moreover, the derivative of p is discontin- 




Figure 4: Eigenvalue density p(x). There are three points where p'(x) is discontinuous, 
associated with supersymmetric M-theory circles. 

uous at the endpoints and at the point x = 0, as predicted by c\, c 2 and c 3 above. The 
Wilson loop calculated from the field theory then agrees with the gravity computation, 
and reads 



log (W) 
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3k 



N 1 ' 2 



(4.22) 



4.3 J\f = 8 super- Yang-Mills with flavour 



In this section we consider a family of M = 2 Chern-Simons-matter theories that 



generalize the mirror to the ABJM theory discussed in section |4.1.2| These were 
discussed in section 4 of [33], having been first introduced in |13j . 

One begins with M = 8 super- Yang-Mills with gauge group U(N), which is the 
theory on iV D2-branes in flat space. In M = 2 language we have three adjoint chiral 
matter fields Xi,X 2 ,X 3 , together with the cubic superpotential TrX 3 [Xi,X 2 ]. To 
this we add matter fields in the fundamental and anti-fundamental representations, 
which breaks the supersymmetry generically to M = 2. More precisely, we add n\ 
fields (q]j ,(jj ), n 2 fields (gj , q j ) and n 3 fields (gj ,gj ), together with the cubic 
superpotential 



W 



Tr 



j=i j=i j=i 



(4.23) 
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so that the theory in section 4.1.2 is simply n\ = 1, n 2 = n% = 0. As for the mirror 
to the ABJM theory, the additional (q, q) fields will arise in type IIA from strings 
stretching from the D2-branes to D6-brane loci. 

In [13] it was shown that the (quantum corrected) moduli space of vacua of these 
theories, for N = 1, may be parametrized by the three coordinates A l5 X 2 , X 3 , together 
with the monopole operators T, T, which satisfy the constraint 

Tf = X^X^X^ . (4.24) 

This defines a Calabi-Yau cone C(Y 7 ) as a hypersurface singularity in C 5 . The M-theory 
circle is straightforward to identify in this case, since by definition the monopole op- 
erators T, T have charges ±1, respectively, under U(1)m, while the Xi are uncharged. 
Notice this implies that the quotient C(Yf)/C* M = C 3 by the complexified M-theory 
circle (defined as a GIT quotient) is simply C 3 , parametrized by A 1 ,A 2 ,A 3 , which 
implies that the type IIA description involves N D2-branes in flat space. Moreover, 
U(1) M Axes T = f = 0, which defines the surface {X^X^X™ 3 = 0} C C 3 , which be- 
comes a D6-brane locus. This then geometrically engineers the gauge theory described 



above, with superpotential (4.23). 

It is straightforward to analyse the matrix model for this gauge theory, as described 
in section 



2.2 and carried out in p3j. The eigenvalue density is given by 



(E, = l n » A i- 2A m) / _ N ^ ^ n 

n ( T \ - ) 87r 2 AiA 2 A 3 ^ Xmm ) ' Xmln ^ X ^ U (A OK) 

87r 2 AiA2A3 l-^max %) > U < X < X mSLX 



and the endpoints are 



, / 8^A 1 A 2 A 3 (Et 1 nA T 2A m ) 
max/mm " V (ELi nA^EU **i ± 2A m ) ' { b) 

Here Aj = A(Xj), i = 1,2,3, are the R-charges of the fields X iy while A m = A(T) = 



A(T) is the R-charge of the monopole operators. As described in section 2.2, these may 
be left a priori arbitrary at this point, the only restriction being that the superpotential 
W has R-charge A(W) = 2. This leads to the constraint ^2 i=1 Aj = 2. The shape of 
p as a function of x is shown in Figure [5] 

The superconformal R-charges are determined by maximizing the free energy F as 
a function of the R-charges. This immediately leads to A m = 0, and then 



2^271 A /A 1 A 2 A3 (V 3 , mAA 

V " ^ ,=1 V / 2 , (4.27) 



3 

\3 



which must be further maximized subject to the constraint Yli=i ^ = 2- In practice 
the formulae are rather too unwieldy for general n,, so following [33] we restrict to the 
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Figure 5: Eigenvalue density as a function of x. There are three points where p'(x) is 
discontinuous, and we correspondingly expect to find three critical points of h M , with 
associated supersymmetric circles. 



case ni = a,ri2 = 6, n^ = 0. In this case the free energy is maximized by 



Ax = ■ 
and thus 



- 26 + Va 2 + b 2 -ab 
2(a - b) 



A 9 



b - 2a + Va 2 + b 2 - ab 
2(6 -a) 



A, 
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2 



x 



max/min 



±2tt 



AiA 2 



aAi + 6A 2 



(4.28) 



(4.29) 



The moduli space equation (4.24) correspondingly reduces to TT = XfX^- The field 
X 3 is then unconstrained, and the Calabi-Yau cone takes the product form C(Y 7 ) = 
C x C(Y 5 ), where X 3 is a coordinate on C and C(Y 5 ) is precisely the Y 5 
singularity^] The toric diagram has lattice vectors 



" a,b,a 



toric 



Wl = (0,0,0), w 2 = (0,1,0), w 3 = (1,0,0) 
w 4 = (0,0, a) , w 5 = (0,1,6) , 



(4.30) 



and is shown in Figure [6j Recall that we parametrize the Reeb vector by £ = 
(4, £2, £3, £4), and that the four- dimensional outward-pointing vectors to the facets are 
v a = (1, w a ). With the method described earlier in section 3.5, the volume of the base 



I7 and the Reeb vector can be found and expressed in terms of Ai and A 2 , and one 
finds 

^ 1 



Voi(r 7 



7T 



6 A 1 A 2 (aA 1 + 6A 2 



and 



£ = (4,l,2A 2 ,aA 1 + 6A s 



(4.31) 



(4.32) 



The M-theory circle in this basis is given by Cm = (0, 0, 0, —1); one can derive this 
by writing the functions T, T, Xi in terms of the toric geometry formalism above (see, 



10 A discussion of the relation between the hypersurface and toric geometry descriptions of this 
three- fold singularity may be found in |45j . 
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w, 



w 5 = (0,l,b) 




(0,1,0) 



Figure 6: Toric diagram corresponding to C(Y 7 ) = {TT = X%X%} x C The apex 
is not an isolated singularity, as one sees from the non-triangular face with vertices 
(1,2,4,5). 

for example, section 4.3 of [15]). Recall also that in this formalism the M-theory 
Hamiltonian function is given by (3.33). Thus in this case we have simply h>M = — 1y^. 
The critical points of Hm must always lie on the boundary of the polyhedral cone, 
which are coordinate singularities, and thus it is easiest to determine this critical set 
using the method described at the end of section |3.5[ We denote the face of the toric 



diagram which has vertices {v a , Vb, v c , . . .} by (a, b, c, . . .). Equation (3.34) then has 
two types of solution: 



\hM\ 







1 



on 



on 



aAi + bA 2 
and correspondingly one has the critical values Hm 

1 



c 3 



-ci 



(2,3,5), (1,3,4), (1,2,4,5) 
(1,2,3), (3,4,5), 

= ^ given by 

and c 2 = . 



(4.33) 



(4.34) 



aA 1 + 6A 2 
Notice here that the face (1, 2, 4, 5) (being non-triangular) corresponds to the S 1 locus 



of L a,b,a conical singularities in Yj. Using the general formula (1.6) we then find that 



these values of Cj precisely match the corresponding positions si,a;2, X3 at which the 



derivative of the eigenvalue density p'(x) is discontinuous. Finally, using (1.4) the 
Wilson loop is 



log (W) 



gravity 



2ty 



AiA 2 



0A1 + &A 2 



-N 1 ' 2 



X T j 



^N 1 ' 2 = log (W) QFT , (4.35) 



where we used (4.29). 
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4.4 £ a2aa Chern-Simons-quivers 

In this section and the next we study two families of examples whose matrix models 
were first analyzed in [37] . 

The M = 2 field theories begin life as low-energy theories on N D2-branes at an L a,b,a 
Calabi-Yau three-fold singularity. This may be simply described as the hypersurface 
{wz = u a v b } C C 4 , where (w,z,u,v) are the coordinates on C 4 . This geometry also 
appeared in the previous subsection of course, but there the M-theory Calabi-Yau four- 
fold was a product C x C(L a,b,a ), whereas here instead C(L a,b,a ) arises as the type IIA 
spacetime. The low-energy theory on the N D2-branes is known from J57J EHl EH], and 
is described by a U(N) a+b quiver gauge theory, with a superpotential W consisting of 
both cubic and quartic terms in the bifundamental and adjoint chiral matter fields. 
Without loss of generality we may take b > a, in which case there are b — a adjoint 
chiral superfields assocated to b — a of the a + b U(N) gauge group factors, and a total 
of 2(a + b) bifundamental fields. We refer the reader to the above references for further 
details of these quiver gauge theories. 

The D2-brane theories become M2-brane theories at a Calabi-Yau four-fold by turn- 
ing on RR flux in the type IIA background, following [T2J and in particular the construc- 
tion in [BJ. Geometrically this fibres the M-theory circle over the base C(L a,b,a ), and 
in the field theory introduces Chern-Simons couplings for the gauge group, described 
by a vector of Chern-Simons levels k = (fci, . . . , k a+ b) = (ki, . . . , kb- a \\kb- a +i, ■ ■ ■ , k a +b), 
where the double bar separates the copies of U(N) with adjoint fields from those with- 
out. This construction is described in more detail in [37J. 

Our first class of examples arise from L a,2a ' a quiver theories, where the vector of 
Chern-Simons levels is k = (0, ... ,0, —2k\\k, k, —k, k, —k, . . . , k, —k, k), with k G Z. 
These theories generalize the model first studied in [21] . The matrix model may be 



solved using the general large N saddle point method described in section 2.2, and one 
finds [37] the eigenvalue density 



Aak-nx(l— A)+fi 



< X < 



16a7r 3 (l-A)A 2 ' 4afc7r(l-A) ^ ^ ^ 2afcvr(2-A) 

P\ X ) :: \ 16a7r 3 A(2-A)(l-A) ' _ 2afc7r(2-A) < X < 2afc7r(2-A) ' (4.36) 

4aknx(l-A)-fi M <r t ^ M 

16a7r 3 (l-A)A 2 ' 2afc7r(2-A) ^ X ^ 4afc7r(l-A) 

where we have defined_3 



2 fcA(2-3A + AT 
/i = Savr^/ I -^ K . (4.37) 

Here the single R-charge variable A parametrizes the R-charges of all the chiral matter 
fields, as in [37]. The eigenvalue density p(x) is shown in Figure [7j 



11 The variable [i arises as a Lagrange multiplier, enforcing that p(x) is a density satisfying (2.131 
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X\ X 



^VQ J\>A "^Tl 



Figure 7: Eigenvalue density as a function of x. There are 4 points x%, x 2 , £3, £4 where 
p'(x) is discontinuous, corresponding to critical points of h M . 



The free energy, as a function of A, is given by 



Son / fcA(l-A)g(2-A)g 3/2 
3 V (4-3A) 



(4.38) 



One may then maximize F to determine the superconformal A, finding the cubic 
irrational 



19- 
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-jp - (431 - lSv 7 ^ 



1/3 



~ 0.319 . (4.39) 



(431-18^417) 

This agrees with the value computed in [ST] , which was for the particular case a — 1. 
Turning to the dual geometry, the Calabi-Yau four-fold that arises as the Abelian 
N = 1 moduli space of these theories has toric data (for k = 1) 



Wl = (0,2a,0), w 2 = (-1,0,0), w 3 = (-1,0,0) 
u>4 = (0, a, a) , W5 = (0, a, — a) , wq = (0,0,0) 



(4.40) 



and with toric diagram shown in Figure [8j The volume of Y 7 may be computed as 



described in section 3.5, and one obtains 



Vol(Y 7 ) 



vr 4 (4 - 3A) 



96a 2 £;A(A - 1) 2 (A - 2) 



with corresponding Reeb vector field 



£ = (4,-4A,2a(2-A),0). 



(4.41) 



(4.42) 



The M-theory circle for this field theory is given in this basis by (m — (0,0,0, 1), 
so that Y 7 is given by a Z^ quotient of the geometry appearing in Figure [HJ with 
Zj. C U(1)m- We may again compute the critical points of the M-theory Hamiltonian 



h M = 2|/ 4 using the method at the end of section 3.5 Equation (3.34) has solutions 
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w 4 =(0,a,a) 



w 6 = (0,0,0) 




w 2 =(-l,a,0) 

w x = (0,2a,0) 



w 5 = (0,a,-a) 
Figure 8: Toric diagram of the L a ' 2a,a Chern- Simons-quiver theories with k — 1. 

associated to the following faces of the toric diagram: 



\hhd\ 
\h M \ 



4a(l-A) 
1 



on (1,4,5,6) 

on (2,3,4), (2,3,5), 

(1,2,4), (1,2,5), (3,4,6), (3,5,6), (4.43) 



on 



2a(2 - A) 
and correspondingly one has the critical values Hm — c % given by 

1 , 1 



Cl 



"Ci 



and c 3 



-c 2 



(4.44) 



4afc(l- A) ' ' " z 2ka{2- A) ' 

Note here that the face (1,4,5,6) describes a singular S 1 locus in Y 7 , and thus al- 
though Km — here, formally leading to zero-action M2-branes, the tangent space is 



singular. Using the general formula (1.6) we then find that these values of q precisely 



match the corresponding positions xi,X2,X3,X4 at which the derivative of the eigen- 
value density p'{x) is discontinuous. Explicitly, the actions of M2-branes wrapped on 
the corresponding calibrated S 1 C Y-j are then 



log (W) 



SM2(C2) 
-Sm2 (C4) 



47r(l-A) VJfc|] _ :<A| 

27r(2-A) v / fei; _ : , Ai 



.A ' N X/2 



(4.45) 



with the latter determining the Wilson loop VEV, and showing that the field theory 
and gravity computations of it agree. 
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4.5 £ a Aa Chern-Simons-quivers 



Our second family within this class are the L a ' b ' a Chern-Simons-quiver theories, with 
the vector of Chern-Simons levels now given by k = (0, . . . , k, — 2k\\k, 0, . . . , 0). One 
finds the eigenvalue density [37] 



4fc7rx(l-A)+^ 



p(x) 



167r 3 (l-A)A((6-2)(l-A)+aA) ' 

M 

16vr 3 (l-A)A(6(l-A)+aA) ' 

4fe7rx(l— A)— /l 

~~ 167r 3 (l-A)A((6-2)(l-A)+aA) 

where we have defined 



4fc7r(l-A) 



J!_ 



2k-rr(b(l-A)+aA) 

M 

2fc7r(6(l-A)+aA) 



< X < 
< X < 



2fc7r(6(l-A)+aA) 



2fc7r(6(l-A)+aA) 



(4.46) 



< x < 



/' 



4/cvr(l-A) 



/' 



8tt 2 



'fcA(l- A) 2 (6(l- A) + aA) 2 
(6-2)(l- A) + aA 



(4.47) 



Again, the R-charge variable A parametrizes the R-charges of all the chiral matter 
fields, as detailed in J37J- The eigenvalue density p(x) is shown in Figure [9j 

p(x)* 




•*1 -^min 



X3 X4 X B 



Figure 9: Eigenvalue density as a function of x. There are again 4 points xi,x 2 ,x 3 ,x 4: 
where p'(x) is discontinuous, corresponding to critical points of h M . 



The free energy, as a function of A, is given by 



Svr / fc(l-A)3A(fe(l-A) + aA) 2 
3 V ((6 + 2)(l-A) + aA) 



(4.48) 



One may then maximize F to find an expression (not presented) for the superconformal 
A that depends on a and b. 

The corresponding Calabi-Yau four-fold that arises as the Abelian N = 1 moduli 
space of these theories has toric data (for k — 1) 

Wl = (0,0,0), w 2 = (1,-1,0), w 3 = (1,1,0), w 4 = (6-1,-1,0), 

w 5 = (6-1,1,0) , w 6 = (6,0,0) , w 7 = (0,0,1) , w 8 = (a, 0,1) , (4.49) 



and with toric diagram shown in Figure 10 The volume of Y 7 may be computed as 



described in section |3.5[ and one obtains 

vr 4 ((6 + 2)(l-A) + aA) 



Vol(F 7 



96A;A(1 - A) 2 (6(l - A) + aA) 2 ' 
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(4.50) 



w 7 = (0,0,1) 



w 2 = (1,-1,0) 



w 4 = (b- 1,-1,0) 




w 8 = (a,0,l) 



w 3 = (1,1,0) 



w 5 = (b-l,l,0) 



>--' w 6 = (b,0,0) 
Figure 10: Toric diagram of the L a,b,a Chern-Simons-quiver theories with k = 1. 

with corresponding Reeb vector 

1 = (4,2(6(1- A) + aA),0,4A) . (4.51) 

The M-theory circle for this field theory is given in this basis by (m — (0, 0, 1, 0), so 



that again Yj is given by a Z^. quotient of the geometry appearing in Figure 10 , with 
Zfe C U(1)m- The M-theory Hamiltonian is Km — 2?/3, and its critical points may be 



computed from equation (3.34), which has solutions on the following faces of the toric 
diagram: 

on (1,2,3,4,5,6),) 

on (2,4,7,8), (3,5,7,8), (4.52) 



Km 



4(1 -A) 



\h 



1 



M\ 



on 



(1,2,7), (1,3,7), (4,6,8), (5,6,8) 



2(6(1- A) + aA) 
and correspondingly one has critical values hu = Q given by 

1 . 1 



Cl 



-Cl 



Ak(l-A) 



and c 3 = — c 2 



2k(b(l - A) + oA) 



(4.53) 



Using the general formula (1.6) we then find that these values of C; precisely match the 



corresponding positions xi,X2,X3,X4 at which the derivative of the eigenvalue density 
p'(x) is discontinuous. Explicitly, the actions of M2-branes wrapped on the correspond- 
ing calibrated S 1 C Y 7 are 

- S M2 (c 2 ) = 4tt(1 - A)J k(l2+h)l t A) ^ N^ , 

log(W) = -S M 2(c A ) : 2;ri/;CI - A ■■■ r/A) ) t / ,.,.,,. ,,, A ,, .." a^ 1 ^ 



(4.54) 
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with the latter determining the Wilson loop VEV, and showing that the field theory 
and gravity computations of it agree. 

4.6 Duals to non-Einstein solutions with flux 

In this final class of examples we examine a family of M = 2 superconformal theories 
which are dual to warped non- Einstein solutions, with non-zero internal _F 4 flux in 



(3.1). As we described in section 3.1. solutions with non-zero M2-brane charge N still 
have a contact structure, and our formalism then applies. 

Following [33], we begin with the M = 8 super- Yang-Mills theories with flavour in 



section 



4J for which the dual Calabi-Yau four-fold is C(Y 7 ) = C x C{L a ^ a ). The 
field theories have a free chiral field X 3 , of scaling dimension A(X 3 ) = A 3 = |, and we 
consider perturbing the theory by adding the deformation AXf to the superpotential. 
In three dimensions this is a relevant deformation for p = 2 and p — 3. The gravity 
dual to the resulting infrared (IR) fixed point of the massive p = 2 deformation is the 
Corrado-Pilch-Waner solution of [20], while the supergravity dual to the cubic p = 3 
deformation was found only recently (independently) in [251 160] . In both cases these 
are warped AdS 4 x Y 7 solutions with flux, where crucially Y 7 has the same topology as 
before the deformation, and the toric symplectic structure of the cone C(Y 7 ) is also the 
same. However, the Reeb vector field (and hence also contact structure) is different. 

This deformation was studied in the matrix model in [33], where they found the 
universal behaviour 

Q - 16( V )3/2 , (4.55) 

relating the free energies F of the IR and UV theories. One correspondingly finds that 
the scaling dimensions A, = A(Xj) are related as 

A} R 4(p-l) A IR aTR 2 tl N 

and that the eigenvalue distribution takes the same form as before, with x max = — x m i n , 
but now with the endpoints rescaled as 



-^max V 3p 



As explained in general in [25], using (2.16) the field theory result (4.55) matches 



with the supergravity result relating the (contact) volumes of the IR and UV solutions: 

VoL,(r 7 IR ) 27p 4 



Vol(Y 7 uv ) ' 256(p - l) 3 
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V(p) . (4.58) 



This result was derived in [23] using the contact geometry of the IR solutions. It is 
crucial that here the volume is the contact volume Vol^ and not the Riemannian volume 
(which is different). The Reeb vector field of the UV Calabi-Yau geometry C x C{Y§) 
may be written as 



-uv 



<9 03 + d v , (4.59) 



where <9^ 3 rotates the free field X3 with charge 1 (so X3 = |X3|e 1(?i3 ), while d v denotes 
the Reeb vector field of the Calabi-Yau three-fold cone C(Y 5 ). As shown in [25], the 
IR Reeb vector field for the solution with flux is then 



.ir 4 Q , 4(p-l) 



e R = -fy,+ \ ' d^. (4.60) 

p 3p 



In fact these formulae are directly related to the rescalings of the R-charges in (4.56). 



Using this simple rescaling one can check that the M2-brane actions are rescaled as 

Sm2 i(FI) 4v^^T 



S$ " y/V{p) " VSp 



(AM) 



thus matching the field theory result (4.57). It of course follows immediately that the 



field theory and gravity results for the Wilson loop VEV agree. What is remarkable 
here is that we are able to compute this so straightforwardly, even though the p = 3 
supergravity solution is known only numerically in [251 160] . Of course, it is precisely be- 
cause we know the contact structure explicitly that the computation is straightforward, 
and this is all that is required to compute these BPS quantities. 

5 Discussion 

In this paper we have shown that the large N field theory and gravity computations of 
the BPS Wilson loop agree in a large class of three-dimensional M = 2 superconformal 
field theories with AdS4 x Y 7 gravity duals. In fact really this matching is a corollary 
of the fact that the image of the M-theory Hamiltonian hMiXi) = [ c min, c max ] is equal 
to the support [x min , x max ] of the real part of the saddle point eigenvalue distribution 
in the large N matrix model, with the proportionality factor between the variables 



x and c given by (1.6). Moreover, the critical points of Km, which give the loci of 
supersymmetric M2-branes wrapping the M-theory circle, always map under Hm to the 
points at which p'(x) is discontinuous in the matrix model. 

In [6T1I62] a relation was also found between p(x), and other matrix model variables, 
and certain geometric invariants. In particular, equation (1.4a) of [62] relates p(x) to 
the derivative of a function that counts operators in the chiral ring of the gauge theory 
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according to their R-charge and monopole charges. In the language of the current 
paper, the monopole charge here is the charge under U(1)m- It would be interesting 
to explore the relation between this formula and the critical points of the Hamiltonian 
function hu- We also note that a similar relation arises when one computes the vacuum 
moduli space of the theory. To see this, recall that the Calabi-Yau four-fold cone 
C(Y?) arises as (a component of) the moduli space of vacua of the Abelian (N = 1) 
field theory, in which <j\ — • ■ ■ — a g = s are all equal. Moreover, the value of s is 
then precisely related to the D-term for the M-theory circle action on C{Yt) - see the 
discussion in section 3 of [3J. This D-term is, in mathematical language, then the M- 
theory Hamiltonian function. We leave further exploration of these issues for future 
work. 
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